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Abstract 

We study the dependence of the heat transfer and the van der 
Waals friction between two semi-infinite solids on the dielectric prop- 
erties of the bodies. We show that the heat transfer and van der 
Waals friction at short separation between the solids may increase by 
many orders of magnitude when the surfaces are covered by adsor- 
bates, or can support low-frequency surface plasmons. In this case 
the heat transfer and van der Waals friction are determined by reso- 
nant photon tunneling between adsorbate vibrational modes, or sur- 
face plasmon modes. The enhancement of the van der Waals friction is 
especially large when in the adsorbed layer there is an acoustic branch 
for the vibrations parallel to the surface like in the case of Cs adsorp- 
tion on Cu(100) surface. In this case we show that even for separation 
d = lOnm, the van der Waals friction induced by adsorbates can be so 
large that it can be measured with the present state-of-art equipment. 
The van an der Waals friction is characterized by a strong distance 
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dependence (~ l/^ 6 ), and at the small distances it can be much larger 
than the electrostatic friction observed in f5]. 

Keywords: non-contact friction, van der Waals friction, radiative 
heat transfer, atomic force microscope, adsorbate vibrational mode 

1 Introduction 

In the history of physics the studies of thermal radiation from materials 
always played a very important role. Here, it is enough to mention that 
quantum mechanics originated from the attempts to explain a paradoxical 
experimental results related to the black body radiation. In the past, the 
nonradiative near-field part of electromagnetic radiation usually was ignored, 
because it plays no role in the far-field properties of emission from planar 
sources. Nevertheless, recent interest in microscale and nanoscale radiative 
heat transfer [U 121 El IH E], together with the development of local-probe 
thermal microscopy jHj have raised new challenges. These topics, and the 
progress in detecting non-contact friction on sub-attonewton level [3 El 13 
[HUE]; an d the observation of coherent thermal emission from doped silicon 
and silicon carbide (SiC) gratings [12] have in common the substantial role 
of the nonradiative (evanescent) thermal field. 

It is well known that the radiative heat transfer between two bodies sep- 
arated by d » At = ch/ksT is given by the Stefan-Boltzmann law: 



where T\ and T2 are the temperatures of solid 1 and 2, respectively. In this 
limiting case the heat transfer between two black bodies is determined by 
the propagating electromagnetic waves radiated by the bodies, and does not 
depend on the separation d. For d < At the heat transfer increases by many 
orders of magnitude due to the evanescent electromagnetic waves that decay 
exponentially into the vacuum; this is often refereed to as photon tunneling. 
At low temperatures (a few K) it is possible for this form of heat transfer to 
be dominant even at spacing of a few mm, see table 1. 

Table 1. Critical distance At as a function of temperature. For surface 
separation d < At the heat transfer is dominated by the contribution from 
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the evanescent electromagnetic modes. At distances of a few nanometers, 
radiative heat flow is almost entirely due to evanescent modes. 

T(K) A r (7WT 

1 2298.8 

4.2 545.2 

100 22.9 

273 8.4 

1000 2.3 

The problem of the radiative heat transfer between two flat surfaces was 
considered some years ago by Polder and Van Hove pQ, Levin and Rytov j2] 
and more recently by Pendry [3| , and by Volokitin and Persson [U E| . In [HI 
HJ E] it was shown that the heat flux can be greatly enhanced if the electrical 
conductivities of the materials are chosen such as to maximize the heat flow 
due to photon tunneling. At room temperature the heat flow is maximal for 
materials with conductivities corresponding to semi-metals. In fact, only a 
thin film (~ 10A) of a high-resistivity material is needed to maximize the 
heat flux [4 . Another enhancement mechanism of the radiative heat transfer 
may arise from resonant photon tunneling between localized states on the 
different surfaces, if the frequency of these modes is sufficiently low to be 
excited by thermal radiation. Recently, it was predicted enhancement of 
the heat transfer due to resonant photon tunneling between surface plasmon 
modes localized on the surfaces of the semiconductors [T31 Ej and adsorbate 
vibrational modes [Bj. 

The problem of radiative heat transfer is closely related to the non-contact 
friction between nanostructures, including, for example, the frictional drag 
force between two-dimensional quantum wells jTJI El Ej , and the friction 
force between an atomic force microscope tip and a substrate jH |HJ El El HI] 

Recently several groups have observed unexpectedly large long-range non- 
contact friction |HJ El E] • The friction force F acting on an atomic 
force microscope tip was found to be proportional to the velocity v , F = Tv. 
At the separation d = 100A the friction coefficient V & 10~ 10 — 10~ 13 kg/s. 
Although the non-contact friction must have an electromagnetic origin, the 
detailed mechanisms is not fully understood yet. 

In a recent Letter, Dorofeyev et. al. [7j claim that the non-contact friction 
observed in |H] is due to Ohmic losses mediated by the fluctuating electro- 
magnetic field . This result is controversial, however, since the van der Waals 
friction has been shown ^7J El to be many orders of magnitude smaller than 
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the friction observed by Dorofeyev et. al. However, we have shown that when 
the surfaces are separated by about lnm, the van der Waals friction can be 
greatly enhanced by resonant photon tunneling ^111 120] ■ Thus, we found 
that resonant photon tunneling between two semiconductors surfaces of SiC, 
which can support surface plasmon modes in mid-infrared region, may give 
rise a 1000-fold (or more) enhancement of the friction, in comparison with 
the case of good conductors. Furthermore, resonant photon tunneling be- 
tween two Cu(100) surface covered by 0.1 monolayer of K may result in six 
orders of magnitude enhancement of friction compared to clean surfaces. At 
the separation lnm we obtained a friction comparable in magnitude with the 
friction observed in the experiment J[|. 

The origin of the van der Waals friction is very closely connected with 
ordinary (attractive) van der Waals interaction. The van der Waals interac- 
tion between atoms (or molecules) arises from quantum fluctuations in the 
electric dipole moment of atoms. The short-lived atomic polarity can induce 
a dipole moment in a neighboring atom or molecule at some distance away. 
The same is true for extended solids, where thermal and quantum fluctuation 
of the current density in one body induces a current density in another body. 
When two bodies are in relative motion, the induced current will lag slightly 
behind the fluctuating current inducing it, this lag is the origin of the van 
der Waals friction. 

In contrast to the van der Waals interaction, for which a well established 
theory exist [21] , the field of van der Waals friction is still controversial. Thus 
different authors [221 [23 12H 123 EH 121] have derived expression for the van 
der Waals friction between two flat surfaces and between a small particle 
and flat surface using different methods, and obtained results, which were 
not confirmed in subsequent calculations [281 U§\ 1301 IT71 152]. 

In ^7j we have developed the theory of the van der Waals friction based 
on a dynamical modification of the well known Lifshitz theory |3T] of the van 
der Waals interaction. In the nonretarded limit, and for zero temperature, 
this theory agree with the results of Pendry |29j . We have also calculated the 
van der Waals friction between two flat surfaces in normal relative motion 
[HH 120] , and found a drastic difference in comparison with parallel relative 
motion. In the limit, when one of the bodies is sufficiently rarefied, our 
theory gives the friction between flat surface and small particle, which agrees 
with the results of Tomassone and Widom 33J . 

The organization of this article is as follows. In Sec. 2 we present a short 
overview of the main principles of the radiative heat transfer, focusing mainly 
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on the adsorbate vibrational mode enhancement of the radiative heat trans- 
fer. In Sec. 3 we study the van der Waals friction, and its enhancement due 
to surface polaritons and adsorbates. In the light of our theoretical results, 
we also discuss non-contact friction experiments. Finally, Sec. 4 contains the 
summary. 



2 Radiative heat transfer 



2.1 Clean surfaces 



According to PjJ EH El the heat transfer between two semi-infinite bodies, 
separated by a vacuum gap with the width d, is given by the formula 

S= / dw(IIi-n 2 )M (2) 
Jo 



where 



M = — dqq 
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+ — / dqqe 2kd x- ^ ' ^ — + [p -> s , 3 

and where the symbol [p — > s] stands for the terms which are obtained from 
the first two terms by replacing the reflection coefficient R p for p-polarized 
electromagnetic waves with the reflection coefficient R s for s- polarized elec- 
tromagnetic waves, and where p — ((cu/c) 2 — g 2 ) 1 / 2 , and k = \p\. The Planck 
function of solid 1 

ni(cj) = hu (e tUj/kBTl -l)' , (4) 

and similar for II2. The contributions to the heat transfer from the propa- 
gating (q < uj/c) and evanescent (q > ui/c) electromagnetic waves are deter- 
mined by the first and the second terms in Eq.fjHJ), respectively. 

From Eq.® it is easy to see that the propagating photon modes give 
the main contribution to the heat transfer for q < A^ 1 , while the evanescent 
modes for q < 1/d. Thus from phase space arguments it follows that the 
number of the channels of the heat transfer available for evanescent waves 
will be by a factor (Xj-/d) 2 larger than the number of the channel available 
for propagating waves. At d = lnm and T = 300K, this ratio is of the order 
~ 10 8 . 
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Let us firstly consider some general consequences of Eq. (J3J). In the case 
of the heat transfer through free photons (q < uj/c), the transfer is maximal 
when both the bodies are perfectly black and have zero reflection amplitude 
R = R r + iRi = 0. The heat flux due to the evanescent waves is a maximal 
when [3] 

Rl + Rl = e 2kd (5) 

Substituting this result into (jSJ) gives the maximal contribution from the 

evanescent waves 

p T 2 2 

fn \evan "'B* He fa\ 

Wzjmax — 24^, ^ ' 

where q c is a cut-off in q, determined by the properties of the material. 
It is clear that the largest possible q c ~ 1/6, where b is an inter- atomic 
distance. Thus, from Eqs.tjHJ) and (0) we get the ratio of the maximum heat 
flux connected with evanescent waves to the heat flux due to black body 
radiation (S z ) max / S BB ~ 0.25 ■ (At/6) 2 . Thus, at room temperature the 
contribution to the heat flux from evanescent waves may be eight orders of 
magnitude larger than the contribution from the black body radiation, and 
the upper boundary for the heat transfer at room temperature: (S z ) max ~ 
10 n Wm- 2 . 

Let us now apply the general theory to a few different materials. For good 
conductors, using Fresnel formulas for reflection coefficient and assuming 
k B T/ATiha < 1, and X T (k B T '/A-ahaf/ 2 < d < \ T (k B T /Ai^ha)- 1 / 2 (a is the 
conductivity) the contribution to the heat transfer from p-polarized waves 
for good conductors (/c#T '/ 'Anha 1) is given by 



s fS o2<MT('MlV /J 

p ' h\ T d \A-nha) 



(7) 



and for the metal with lower conductivity ( k B T / AitTio < 1), for d < Xx{k B T ' j 'Airha)' 
we get 

SpK0 .12^(-^)7l + ln^y (8) 
p hd? \Anha J \ k B T J v ; 

For good conductors the heat flux depends on the separation as ~ 
and increases with decreasing conductivity as cr -1 / 2 . For k B T/Anha < 1 
the heat flux decreases with separation as d~ 2 , and increases with decreas- 
ing conductivity as cr -2 , while the s-wave contribution is distance indepen- 
dent, S s ~ 0.25k B Ta/ X 2 -. Fig. la shows the heat transfer between two 
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semi-infinite silver bodies separated by the distance d, at the temperatures 
Ti = 273K and T 2 = OK. The s- and p-wave contributions are shown sepa- 
rately. The p-wave contribution has been calculated using non-local optics, 
i.e. including spatial dispersion of the dielectric function (the dashed line 
shows the result using local optics). It is remarkable how important the 
s-contribution is even for short distances. Note from Fig. la that the local 
optics contribution to S p depends nearly linearly on \ jd in the studied dis- 
tance interval, and that this contribution is much smaller than the s-wave 
contribution. Both these observations agree with analytical formulas given 
above. However, for the very high-resistivity materials, the p-wave contribu- 
tion becomes much more important, and a crossover to a 1 /^-dependence 
of S p is observed at short separations d. This is illustrated in Fig. lb and lc, 
which show results for the same parameters as in Fig. la, except that the 
electron mean free path has been reduced from / = 560 A (the electron mean 
free path for silver at room temperature) to 20 A (roughly the mean free 
path in lead at room temperature) (Fig. lb) and 3.4 A (of order the lattice 
constant, representing the minimal possible mean free path) (Fig. lc). Note 
that when / decreases, the p-wave contribution to the heat transfer increases 
while the s-wave contribution decreases. Since the mean free path cannot be 
much smaller than the lattice constant, the result in Fig. lc represent the 
largest possible p-wave contribution for normal metals. However, the p-wave 
contribution may be even larger for other materials, e.g., semimetals, with 
lower carrier concentration than in normal metals. For high resistivity mate- 
rials, when hsT/AnTia > 1 the heat flux is proportional to the conductivity 

S p » 0.2^ (9) 

Figure 2 show the thermal flux as a function of the conductivity of the solids. 
Again we have assumed that one body is at zero temperature and the other 
at T = 273K. The surfaces are separated by d — 10 A. The heat flux for 
other separations can be obtained using scaling ~ 1/d 2 which holds for high- 
resistivity materials. The heat flux is maximal at a = 1316(f2- m) -1 . Finally, 
we note that thin high-resistivity coatings can drastically increase the heat 
transfer between two solids. This is illustrated in Fig.3, which shows the 
heat flux for the case when thin films (~ 10 A) of high resistivity material, 
p — 0.14 • cm, are deposited on silver, (a) and (b) shows the p- and s- 
contributions, respectively. Also shown are the heat flux when the two bodies 
are made from silver, and from the high resistivity material. It is interesting 
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to note that while the p-wave contribution to the heat flux for the coated 
surfaces is strongly influenced by the coating, the s-contribution is nearly 
unaffected. 

Another case where the heat transfer can be large is when resonant photon 
tunneling occurs between surface states localized on the different surfaces. 
The resonant condition corresponds to the case when the denominators in 
Eq.© are small. For two identical surfaces and Ri << 1 < R r , where Ri 
and R r are the imaginary and real parts of the reflection coefficient, this 
corresponds to the resonant condition 

R r exp(-qd) = ±1 (10) 

This condition can be fulfilled in spite of the factor exp(— 2qd) < 1 because 
for evanescent electromagnetic waves there is no restriction on the magni- 
tude of real part or the modulus of R. This opens up the possibility of 
resonant denominators for ^> 1. Close to the resonance we can use the 
approximation 

R= ^— , (11) 

UJ — UJq — 11] 

where uo a is a constant. Then from the resonant condition (R r = ±e qd ) we 
get the positions of the resonance 

u± = uj ± io a e~ qd . 

For the resonance condition to be valid the separation Auj = \uj + — uj_\ 
between two resonances must be greater than the width 77 of the resonance. 
From this condition we get that the two poles approximation is valid for 
q < q c ~ ln(2uj a /r))/d. 

For uj > uj a and q c d > 1, we get 

s , ± = ^[n 1 (wb)-n a (w b )]. (12) 

Note, that the explicit d dependence has dropped out of Eq. (JTJJ). How- 
ever, S may still be d- dependent, through the d- dependence of q c . For small 
distances one can expect that q c is determined by the dielectric properties of 
the material, and thus does not depend on d. In this case the heat transfer 
will also be distance independent. 
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Resonant photon tunneling enhancement of the heat transfer is possible 
for two semiconductor surfaces which can support low-frequency surface plas- 
mon modes in the mid-infrared frequency region. The reflection factor R p for 
clean semiconductor surface at a < At |e| is given by Fresnel's formula 

R P = ^T[, (13) 

where e is the bulk dielectric function. As an example, consider two clean 
surfaces of silicon carbide (SiC). The optical properties of this material can 
be described using an oscillator model e 34 

<«>) =^{l+ / L ~/ T . v ) (14) 

with eoo = 6.7, tu L = 1.83 • lO 14 ^ 1 , uj t = 1.49 • 10 14 s-\ and T = 8.9 • lO 11 ^ 1 . 
The frequency of surface plasmons is determined by condition e r (u) p ) = — 1 
and from (|14|) we get u p = 1.78 • 10 14 s _1 . The resonance parameters 

uj 2 — uj 2 

u a = -^ £ = 8.2 ■ 10 12 s _1 , rj = T/2, q c = 3.6/ d, and cu w u p 

Using these parameters in Eq. (jl2j) and assuming that one surface is at 
temperature T = 300 K and the other at K we get the heat flux S(d) 
between two clean surfaces of silicon carbide (SiC): 

S p w 8.4 ■ 10 9 ^W- m" 2 (15) 

where the distance d is in Angstr0m. Note that the heat flux between the 
two semiconductor surfaces is several order of magnitude larger than between 
two clean good conductor surfaces (see Fig.l). 

2.1.1 Adsorbate vibrational mode enhancement of the radiative 
heat transfer 

Another resonant photon tunneling enhancement of the heat transfer is pos- 
sible between adsorbate vibrational modes localized on different surfaces. Let 
us consider two particles or adsorbates with the dipole polarizabilities ot\(uS) 
and otziu)) and with the fluctuating dipole moments p{ and p 2 normal to the 
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surfaces. Accordingly to the fluctuation-dissipation theorem [33], the power 
spectral density for the fluctuating dipole moment is given by 



n n 

where the Bose-Einstein factor 



(PiPj)u, = ^ Q + ni(u)^j haoti(w)8ij (16) 



nAu) = = ,, „ . (17) 

*v ' e hLu/k B Ti _ ]_ V 7 

Assume that the particles are situated opposite to each other on two different 
surfaces, at the temperatures T\ and T 2 , respectively, and separated by the 
distance d. The fluctuating electric field of a particle 1 does work on a 
particle 2. The rate of working is determined by 

/•CO 

P 12 = 2 dwulma 2 (w)(EuE 12 ) u (18) 
Jo 

where E\ 2 is the electric field created by a particle 1 at the position of a 
particle 2: 

F 8pfl/d3 
12 1 - ai a 2 (8/d 3 ) 2 { ' 

From Eqs. (|16H18j) we get Pi 2 , and the rate of cooling of a particle 2 can be 
obtained using the same formula by reciprocity. Thus the total heat exchange 
power between the particles is given by 

2h r°° , Imailma 2 (8/(i 3 ) 2 . . . . . . . . 

P = Pl2 _P 21 = _ duu ^ Th^-naH) 20 

7r Jo \1- (8/rf 3 ) 2 aia 2 r 

Let us firstly consider some general consequences of Eq. (|2T)|) . There 
are no constrain on the particle polarizability a(u) = a' + ia" other than 
that a" is positive, and a' and a" are connected by the Kramers-Kronig 
relation. Therefore, assuming identical surfaces, we are free to maximize the 
photon-tunneling transmission coefficient 



t= (2i) 



ia"/d 3 y 
\l-(8a/d 3 y 

This function has a maximum when 

a' 2 + a" 2 = (d 3 /8) 2 (22) 
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so that t = 1/4. Substituting this result in (|2U|l gives the upper bound for 
the heat transfer power between the two particles or adsorbates 



Kk 2 2 2 

Pmax — (^l — T 2 ) (23) 

For adsorbed molecules at the concentration n a = 10 19 m -2 , when one surface 
is at zero temperatures and the other is at the room temperature, the max- 
imal heat flux due to the adsorbates S max = n a P max = 10 12 Wm~ 2 , which is 
nearly 10 order of magnitude larger than the heat flux due to the black body 
radiation, Sbb = &bT = 4 ■ 10 2 Wm -2 , where <jb is the Boltzmann constant. 
The conditions for resonant photon tunneling are determined by equation 

a'(oj±) = ±d 3 /8 (24) 

Close to resonance we can use the approximation 

a « , (25) 

UJ — LUq — IT] 

where c = e* 2 /2Muj , and where e* and M are the dynamical charge and 
mass of the adsorbate, respectively. Then from the resonant condition (|2*4*|) 

we get 

u± = ujq ± 8c/ d 3 . 

The separation between the resonances, Au = \u + — U-\ must be greater 
than the width rj of the resonance, so that 8c/ d 3 > rj. 
For rj « 8c/ d 3 , from Eq. (|2UJl we get 

P = ^Mm(u; + ) - n 2 (u + )) + (+ - -)]., (26) 

Using Eq. we can estimate the heat flux between identical surfaces 

covered by adsorbates with concentration n a : S ~ n a P. Interestingly, the 
explicit d dependence has dropped out of Eq. (|2fij) . However, P may still 
be d- dependent, through the d- dependence of u±. For hu± < ksT the 
heat transfer will be only weakly distance dependent. For 8c/c/ 3 < rj we 
can neglect multiple scattering of the photons between the particles, so that 
the denominator in the integrand in Eq. (|2Uj) can be approximated with 
unity. For d » b, where b is the interparticle spacing, the heat flux between 
two surfaces covered by adsorbates with concentration n a \ and n a2 can be 
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obtained after integration of the heat flux between two separated particles. 
We get 

g __ ai a2 / ^ u}lm.ailma2[ni(u) — n 2 (o;)] (27) 

a 4 Jo 

Assuming that a can be approximated by Eq. (|23|) . for u « r\ Eq. (|2Tj) 
gives the heat flux between two identical surfaces: 

^ 'y ^M-^Mi (28) 

For the K/Cu(001) system w = 1.9 ■ 10 13 s _1 , and at low coverage e* = 
0.88e which gives c = e* 2 /2Mu = 7 ■ 10 -17 mV 1 . For r? = lO^s -1 , 
when one surface has T = 300K and the other T = OK, for d > b and 

8a/ d 3 < r] we get 

77 2 

5^5.6-10" 24 -fW- m- 2 (29) 

Cfc 

where the d is in Angstr0m. 

We note that Eq. (|2*H|) can be obtained directly from the heat flux between 
two semi-infinite solids determined by Eq.Q, since in the limit d > b we 
can use a macroscopic approach, where all the information about the optical 
properties of the surface is included in reflection coefficient [6 . The reflection 
coefficient R p , which take into account the contribution from an adsorbate 
layer is given by |¥Tj : 

_ 1 - s/qe + 47m a [sa\\/e + qa±} - qa(l - ATm a qa\\) 
1 + s/qe + 47T77 a |sa||/e — qa±\ + qa(l + 47m a ga||) 

where s = \j q 2 — (w/c) 2 e, and a\\ and a±_ are the polarizabilities of adsorbates 
in a direction parallel and normal to the surface, respectively, e = 1 + 
4nia/u! is the bulk dielectric function, where a is a conductivity, and 77 a is 
the concentration of adsorbates. Eq. (J3Uj) takes into account that the centers 
of the adsorbates are located at distance a away from image plane of the 
metal. Although this gives corrections of the order qa 1 to the reflection 
amplitude, for parallel adsorbate vibrations on the good conductors (when 
e ^> 1), in some cases they give the most important contribution to the energy 
dissipation. As illustration of this macroscoscopic approach, in Fig. (4) we 
show the p-wave contribution to the heat flux for the two Cu(100) surfaces 
covered by a low concentration of potassium atoms (n a = 10 18 m -2 ) and the 
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two clean Cu(100) surfaces. At separation d = lnm the heat flux between 
two surfaces covered by adsorbates is enhanced by five and three orders of 
magnitude in comparison with the p- and s- wave contributions to the heat 
flux between clean surfaces, respectively, and by seven orders of magnitude 
in comparison with the blackbody radiation. 

For d < b the macroscopic approach is not valid any more and we must 
sum the heat flux between each pair of the adatoms. For rj = 10 12 s _1 and 
d < 10A, when one surface has T = 300K and the other T = OK, from 
Eq.()26J) we get the distance independent P ~ 10 _9 W . In this case, for n a = 
10 18 m -3 the heat flux S ~ Pn a rs 10 9 Wm" 2 . Under the same conditions 
the s-wave contribution to the heat flux between two clean surfaces S c i ean ~ 
10 6 Wm~ 2 . Thus the photon tunneling between the adsorbate vibrational 
states can strongly enhance the radiative heat transfer between the surfaces. 
However this enhancement of the heat flux disappears if only one of the 
surfaces is covered by adsorbates. 

It is interesting to note that in the strong coupling case (8c/ d 3 ^> 77) the 
heat flux between two molecules does not depend on the dynamical dipole 
moments of the molecules (see Eq. (J26|) . However, in the opposite case of the 
weak coupling (8c/ d 3 <C 77) the heat flux is proportional to the product of 
the squares of the dynamical dipole moments (see Eq. (j2*5|) ). 



3 Van der Waals friction 

3.1 Clean surfaces 

The frictional stress <r±(\\) which act on the surfaces of two bodies in normal 
(parallel) relative motion can, to linear order in sliding velocity v , be writ- 
ten in the form: cr_i_(ii) = r Y±n\)V. For bodies in parallel relative motion at 
separation d << the friction coefficient 711 is given by ^7] 

x lmR lp lmR 2p — — -— ^ + [p - s}. (31) 
When the two bodies move toward or away from each other the friction 
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coefficient is given by [THJ |20] 



1X Z JO \ OUO J Ju/c 

x [(lmR lp + e~ 2kd \R lp \ 2 lmR 2p )(lmR 2p + e~ 2kd \R 2p \ 2 lmR lp ) 

+ e- 2M (Im(iZ lp i^)) 2 ] _ 2fc ^ ^ 4 + [p - s], (32) 
where the Bose-Einstein factor 



e huj/k B T _ I ' 



At resonance the integrand in Eqs. (|32|) has a large factor ~ ^/R 2 , in 
sharp contrast to the case of parallel relative motion (see Eq.([3*T|)) , where 
there is no such enhancement factor. Thus, at resonance if R 2 << 1 the 
friction for normal relative motion will be much larger than the friction for 
parallel relative motion. In contrast to the heat transfer, the van der Waals 
friction is very sensitive to presence of low frequency excitations which absorb 
plenty momentum without absorbing much energy. Thus the van der Waals 
friction is very sensitive to the type of the material. 

Assuming that the medium 2 is sufficiently rarefied and consist of particles 
with the radius r « d, with the polarizability a (a;) given by Eq.([3*T|). it is 
easy to calculate friction coefficient between a small particle and a flat surface: 



- f°°duu(~) r dqq 2 e- 2 ^ d lma 2 
IX Jo \ OLU J Jo 

x < 2Imi?i p 




+ ^ lmR ls (u) } (33) 

and = 2ry. In the nonretarded limit (c — * oo) this formula agrees with 
the result obtained in [3*3*] . 

For good metals {k B T/A-nha » 1) using for \ T {k B T /A-nhaf' 2 < 
d < \T{^ha/{k B Tfl 2 (A T = ch/{k B T)), we get 

h ( k B T\ 1/2 , , 
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and for d < \ T (k B T / 4:iiha) 3/2 we get 



2 




(35) 

The last contribution will dominate for metal with not too high conductivity 
(k B T/4nha ~ 1). 

For comparison, the p-wave contribution for parallel relative motion for 
d < A c , (A c = c/^vra^T) 1 / 2 ) is given by [T7] 

It is interesting to note that, in contrast to parallel relative motion, for 
normal relative motion of good conductors, for practically all d > the main 
contribution to friction comes from retardation effects, since Eq. f!34|) . in 
contrast to Eq. (J36|) . contains the light velocity. 

From Eq. (|3*2*|) we get the s-wave contribution to friction for d < A c 

7IT" ^10- 2 ^(3-51n(2d/A c )) (37) 

For parallel relative motion the s-wave contribution is two times smaller:7^ an - 

r\ evan 
Z l\\s ■ 

Figures 5 and 6 show the calculated contribution to the friction coef- 
ficient 7 from the evanescent electromagnetic waves for two semi-infinite 
solids, with parameters chosen to correspond to copper (r _1 = 2.5 ■ 10 13 s _1 , 
Up = 1.6- lO^s^ 1 ) at T = 273 A', for parallel (Fig.5) and normal (Fig.6) rela- 
tive motion. Results are shown separately for both the s- and p- wave contri- 
bution. The dashed line shows the result when the local (long-wavelength) 
dielectric function e(uj) is used, and the full line show the result obtained 
within the non-local optic dielectric formalism, which was proposed some 
years ago for the investigation of the optical properties of metals in anoma- 
lous skin effect frequency region |57j . Fig.5 shows that, for sufficiently small 
separations (d < 1000 A), for parallel relative motion the non-local optic ef- 
fects become important for the p— wave contribution. However, for the s— 
wave contribution, for both parallel and normal relative motion, the non- 
local optic effects are negligibly small for practically all the separations. For 
normal relative motion, for the p-wave contribution the non-local optic effects 
are less important, than for the parallel relative motion. 
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For high-resistivity metals (kBT/Airha > 1) for d < A c we get 

lx ~ 0.48——^— 38 
' d 4 Anker K ' 

and 7|| ~ 0.l7 perp . Thus, in contrast to the heat flux, the van der Waals 
friction diverges in the limit a — > 0. Of course, in reality the friction must 
vanish in this limit since the conductivity is proportional to concentration 
of free electrons, and the friction must vanish as the carrier concentration 
vanishes. The origin of the discrepancy lies in the breakdown of the macro- 
scopic theory which was used in the calculation of friction at low electron 
concentration. The macroscopic approach for the electromagnetic properties 
of material is valid only when the length scale of the spatial variation of 
the electromagnetic field is much larger than the average distance between 
the electrons. For evanescent waves this length scale is determined by the 
separation d between the bodies. Thus, the macroscopic approach is valid if 
d >> 72 -1 / 3 , where n is the concentration of electrons. This fact was over- 
looked in Ref . |4Uj . From this requirement we can estimate the maximum 
friction which can be obtained for high resistivity metals. The minimum 
conductivity can be estimated as 



e 2 r 



a 



and the maximum of friction 



d 3 m 



d A AnTio m - m Ane 2 rd 



To estimate the friction coefficient V for an atomic force microscope tip 
with radius of curvature r >> d we can use an approximate formula 



POO 

2tt / d PP1 (z(p)) (39) 
Jo 



10 

where it is assumed that the tip has cylinder symmetry. Here z(p) denotes the 
tip - surface distance as a function of the distance p from the tip symmetry 
axis, and the friction coefficient j(z(p)) is determined by the expressions 
for flat surfaces. We assume that the tip has a paraboloid shape given [in 
cylindrical coordinates (z,p)] by the formula: z = d + p 2 /2r. If 



7(P)= 7- C ^„ (10) 
(d + 



£ 

'IT- 



IS 



we get 

_, 2nr C 2irrd ... , ... 
n — ld n 1 n — 1 
where A e ff is the effective surface area. For high-resistivity metals n — 4, 
and the maximum friction coefficient for spherical tip: 

mk B Tr 

r max ~ lma X dr ~ — — 

Using this formula for r ~ 10~ 15 s, r ~ l/xm and T = 300K we get r max ~ 
10~ 15 kg/ s. This friction is two order of magnitude smaller than was observed 
in a recent experiment [§] at d — lOnm. In the case of the cylindrical tip 
with the width w: 

^ r r— mk B TwR - 5 

r max ~ ImaWdrW ~ 47re2rrfQ , 5 

For w = 7fim and d = lOnm the friction is of the same order as it was 
observed in experiment. Thus, van der Waals friction between high resistivity 
material can be measured with the present state-of-art equipment. 

As in the case of the radiative heat transfer, the van der Waals friction can 
be greatly enhanced when resonant photon tunneling between localized sur- 
face states, e.g. surface plasmon polaritons and adsorbate vibration modes, 
occurs. Using the same approximation as when deriving Eq. p2jl for normal 
relative motion we get 

7± = -li!^ I (41) 

Similar, for parallel relative motion 

7|I = ^ \ m 

l287rk B Tsmh 2 {huJo/2k B T) v ; 

where q c = min(6, ln(2u a /r))/d), and where b is of the order of an interatomic 
distance. Thus if hi(2u a /r)) / d > b, Eq. (|l2*jl is independent of the distance. 
For small distances one can expect that q c is determined by the dielectric 
properties of the material, and does not depend on d. In this case the friction 
will be also distance independent. 

Resonant photon tunneling enhancement of the van der Waals friction 
is possible for two semiconductor surfaces which can support low-frequency 
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surface plasmon modes. As an example we consider two clean surfaces of 
silicon carbide (SiC). The optical properties of this material were described 
above. Using the same parameters as before and at T = 300K we get j\\ ~ 
(10 3 /(i 4 )kg- s~ 1 m~ 2 , where the distance d is in Angstr0m, and 7^ ~ 37M. 
Note that the friction between the two semiconductors is about three order 
of magnitude larger than between two clean good metallic conductors (see 
Fig.(5,6)). 

3.2 Adsorbate vibrational mode enhancement of the 
van der Waals friction 

Another enhancement mechanism of van der Waals friction is connected with 
resonant photon tunneling between adsorbate vibrational modes localized on 
different surfaces. In [HJJ I2H1 we have shown that resonant photon tunnel- 
ing between two surfaces covered by a low coverage of potassium atoms at 
d = lnm gives rise to enhancement in friction by six orders of the magni- 
tude in the comparison with friction between clean surfaces. The adsorbate 
induced enhancement of van der Waals friction is even larger for the case of 
Cs adsorption on Cu(100). In this case even at small 0.08) in 

adsorbed layer there is acoustic branch for vibrations parallel to the surface 
36J . In this case at small frequencies the reflection coefficient is given by [H] 

2qau 2 , 

UJ Z — UJ Z + IUJT] 

where u 2 = 4irn a e* 2 aq 2 /M, r\ is the damping constant for the adsorbate 
vibrations parallel to the surface, a is the separation between the adsorbate 
center and image plane. Using Eq.()43|) in Eq.([31|) for 



a jAnn a e* 2 a 
^d\ Md 2 < 1 

we get 

711 - °- 62 ^T ( 44 ) 

It is interesting to note that the dependence on n a , e*, and M is dropped 
out from Eq. (|44j) . However, Eq.()43|) is only valid when there are acoustic 
vibrations in the adsorbed layer. For Cs adsorbed on Cu(100) surface the 
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acoustic vibrations exist only for 9 > 0.1 The friction coefficient for the 
atomic force microscope can be estimated using approximate Eq. (J39j) . Using 
Eq.(J33J for a cylindrical tip we get 

r; „ o.6S kBT '*y (45) 

II ^5.5 V / 

where r is the radius of the curvature of the tip and w is its width. In the case 
of Cs adsorbed on Cu(100) surface the damping constant rj ~ 3 ■ 10 9 s _1 and 
a = 2.94A gT]. Than for r = 1/im, w = 7/j.m, T = 293 K at d = lOnm we get 
T|| = 0.5-10 _13 kg/s, that is only three times smaller than the friction observed 
in [5| at the same distance. However van der Waals friction is characterized 
by much stronger distance dependence (~ 1/d 5 ' 5 ) than in experiment (~ 
1/cf 1 , where n = 1.3 ± 0.2). Thus at smaller distances the van der Waals 
friction will be much larger than friction observed in [Hj and can be measured 
experimentally. Fig. 7 shows the friction coefficient between the copper tip 
and the copper substrate as a function of the separation d, when the surfaces 
of the tip and the substrate are covered by low concentration of the Cs 
atoms and for the clean surfaces. In comparison, the friction between two 
clean surfaces at the separation d = lnm is eleven orders of the magnitude 
smaller. However, the friction between clean surfaces shown in Fig. (7) was 
calculated in local optic approximation. For parallel relative motion the non- 
local optic effects are very important (see Fig. 5) and when these non-local 
optic effects are taken into account the friction between adsorbate covered 
surfaces at d — lnm will be by seven orders of magnitude larger than the 
friction between clean surfaces in parallel relative motion. 



4 Summary 

We have studied how the radiative heat transfer and van der Waals friction 
between two bodies depends on the dielectric properties of the media. We 
have found that, at short distances between the bodies, the thermal flux 
can be significantly enhanced in comparison with the black body radiation, 
in particular when the material involved can support low-frequency adsor- 
bate vibrational modes, or surface plasmon modes, or the conductivity of 
the metals is chosen to optimize the heat transfer. This fact can be used in 
the scanning probe microscopy for local heating and modification of surfaces. 
We have shown that the van der Waals friction can be enhanced by several 
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orders of magnitude in the case of resonant photon tunneling between low- 
frequency surface plasmon modes and adsorbate vibrational modes. In the 
case of friction between two Cu(100) surfaces covered by a low concentration 
of cesium atoms at d — lOnm we have found the friction to be of the same 
order of the magnitude as it was observed in experiment [9 . However, the 
van der Waals friction is characterized by stronger distance dependence than 
the friction observed in experiment. Thus at small distances the van der 
Waals friction can be much larger than the friction observed in jS] and can 
be measured experimentally. This effect can be important technologically 
for ultrasensitive force registration and from basic point of view. The fric- 
tion observed in [5] can be explained by electrostatic friction jjl] when the 
electromagnetic field in the vacuum gap is mediated by bias voltage or by 
inhomogeneities of the surfaces. 
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FIGURE CAPTIONS 

Fig. 1 (a) The heat transfer flux between two semi-infinite silver bodies 
as a function of the separation d, one at temperature T\ = 273 K and another 
at T 2 = K. (b) The same as (a) except that we have reduced the Drude 
electron relaxation time r for solid 1 from a value corresponding to a mean 
free path vpT = I = 560 A to 20 A. (c) The same as (a) except that we 
have reduced I to 3.4 A. For silver (Fig. la) at T = 273 K the conductivity 
a = 5.6 • lO^s" 1 and k B T/Arrfia = 4.6 ■ 10" 6 and for Fig. lb and Fig.2c these 
quantities can be obtained using scaling a ~ Z -1 . (The base of the logarithm 
is 10) 

Fig. 2 The thermal flux as a function of the conductivity of the solids. 
The solid surfaces are separated by d = 10 A. The heat flux for other sepa- 
rations can be obtained using scaling ~ 1/d 2 which holds for high- resistivity 
materials. (The base of the logarithm is 10) 

Fig. 3. The heat flux between two semi-infinite silver bodies coated with 
10 A high resistivity (p = 0.14 ficm) material. Also shown is the heat flux 
between two silver bodies, and two high-resistivity bodies. One body is at 
zero temperature and the other at T = 273A. (a) and (b) shows the p and 
s-contributions, respectively. (The base of the logarithm is 10) 
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Fig. 4. The heat flux between two surfaces covered by adsorbates and 
between two clean surface, as a function of the separation d. One body is at 
zero temperature and the other at T = 273 K. For parameters corresponding 
to K/Cu(001) and Cu(OOl) [M] {u± = 1.9 ■ lO 13 ^ 1 ,^ = 4.5 • lO 12 ^ 1 ,^ = 
2.8 ■ 10 10 s _1 , 771 = 1.6 ■ 10 12 s -1 , e* = 0.88e) (The base of the logarithm is 10) 

Fig. 5. The friction coefficient for two flat surfaces in parallel relative 
motion as a function of separation d at T — 273 K with parameter cho- 
sen to correspond to copper (r _1 = 2.5 • 10 13 s _1 , u p = 1.6 • 10 16 s _1 ). The 
contributions from the s— and p— polarized electromagnetic field are shown 
separately. The full curves represent the results obtained within the non- 
local optic dielectric formalism, and the dashed curves represent the result 
obtained within local optic approximation. (The base of the logarithm is 10) 

Fig. 6. The friction coefficient for two flat surfaces in normal relative 
motion as a function of separation d at T — 273 K with parameter cho- 
sen to correspond to copper (r -1 = 2.5 • 10 13 s _1 , u p = 1.6 ■ 10 16 s _1 ). The 
contributions from the s— and p— polarized electromagnetic field are shown 
separately. The full curves represent the results obtained within the non- 
local optic dielectric formalism, and the dashed curves represent the result 
obtained within local optic approximation. (The base of the logarithm is 10) 

Fig. 7. The friction coefficient between the copper tip and copper sub- 
strate which surfaces are covered by low concentration of cesium atoms, as a 
function of the separation d. The cylindrical tip is characterized by radius of 
curvature r = 1/im and width w = 7/im. For other parameters correspond- 
ing to Cs adsorbed on Cu(100) surface at coverage 6 « 0.1 and for Cu(100) 
HHE]: e* = 0.28e, 77 = 3 • 10V 1 , a = 2.94A, T = 293 K. (The base of the 
logarithm is 10) 
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